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Abstract 


Understanding the unsteady aspects of turbine rotor flow fields is critical to successful 
future turbine designs. A technology program was conducted at NASA’s Marshall Space Flight 
Center to increase the understanding of unsteady environments for rocket engine turbines. The 
experimental program involved instrumenting turbine rotor blades with miniature surface- 
mounted high frequency response pressure transducers. The turbine model was then tested to 
measure the unsteady pressures on the rotor blades. The data obtained from the experimental 
program is unique in two respects. First, much more unsteady data was obtained (several 
minutes per set point) than has been possible in the past. Also, an extensive steady performance 
database existed for the turbine model. This allowed an evaluation of the effect of the on-blade 
instrumentation on the turbine’s performance. A three-dimensional unsteady Navier-Stokes 
analysis was also used to blindly predict the unsteady flow field in the turbine at the design 
operating conditions and at +15 degrees relative incidence to the first-stage rotor. The predicted 
time-averaged and unsteady pressure distributions show good agreement with the experimental 
data. This unique data set, the lessons learned for acquiring this type of data, and the 
improvements made to the data analysis and prediction tools are contributing significantly to 
current Space Launch Initiative turbine airflow test and blade surface pressure prediction efforts. 

Nomenclature 

Co: Spouting Velocity (i t/s) 

C p : Specific heat at constant pressure (BTU/lb m R) 

FP: Flow parameter 

g c : Conversion constant (32.174 (ft*lb m )/(lbf*s‘)) 

J: Conversion constant (778.3 ft-lbf/BTU) 
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M: Mach number 

N: Speed (RPM) or synchronous spectral component 

P: Pressure (psia) 

Pr: Pressure ratio 

T: Temperature (°R) 

Tq: Torque (ft-lbf) 

U: Disk tangential speed (ft/s) 

U: Uncertainty 

W : Mass flow rate (lb n /sec) 

a: Absolute circumferential flow angle measured from the axial direction 

P: Relative circumferential flow angle measured from the axial direction 

y: Ratio of specific heats 

q: Efficiency 

Subscripts 

0: Total 

1: Inlet 

2: Exit 

ave: Average 

Fac: Facility 

max: Maximum value 

min: Minimum value 

inst: Instantaneous value 

th: Thermodynamic method 



t-s: total-to-static 

t-t: total-to-total 

Introduction 

The goals of future and next generation reusable launch systems are to increase safety and 
reliability, to reduce unit and/or operational costs (life, time between replacement and overhauls, 
operations complexity), and to reduce weight. To meet these goals, rocket engine components 
are required to be smaller, lighter weight, higher performing, more reliable, and less costly. 
These requirements push turbine designs to operate beyond the limits of the past. The flow 
through a turbine stage has always been extremely complicated. The rotor flow field in particular 
is unsteady and generally ihree-dimensional. This flow field unsteadiness is a major factor in 
turbine performance and life, and, as the turbine designs become more compact and closely 
coupled to meet new size and weight requirements, the flow field unsteadiness increases. 
Therefore, understanding the unsteady aspects of the flow field is critical to the successful design 
of future high performance turbines. The fundamental understanding of the flow physics will 
lead to more efficient designs that could significantly reduce the unsteadiness to provide more 
durable, higher performing turbines to meet program goals. 

A technology program was conducted at Marshall Space Flight Center (MSFC) to better 
understand the unsteady environments of rocket engine turbines. The program involved 
instrumenting the 1 st stage rotor blades of the space shuttle main engine (SSME) high-pressure 
fuel turbopump (HPFTP) turbine with surface-mounted high frequency response pressure 
transducers. The HPFTP turbine model was then tested in air in the MSFC turbine test facility 
(TTE) to measure the unsteady pressures on the rotor blades. In addition, a time-accurate 
computational fluid dynamics (CFD) analysis was performed in a computationally efficient 
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manner to predict the unsteady flow environment using the SSME blade surface pressure data for 
validation. 

This article discusses the experimental program and unsteady CFD results. It describes 
the air test conducted in the TTE including sections describing the test facility, the test model, 
and both the steady and unsteady instrumentation and data acquisition equipment. Issues critical 
to the success of the test, such as transducer calibration and data rates, are discussed. Steady 
state turbine performance data is presented and compared with performance data obtained from a 
previous test of the model in the same configuration with no unsteady instrumentation. Blade 
surface pressure data is then presented to demonstrate data repeatability and validation. Next, 
reduced time averaged and lime resolved pressure data are discussed. Finally, the unsteady CFD 
methodology and results aie presented along with comparisons to the experimental data. This 
unique data set, the lessons learned for acquiring this type of data, and the improvements made to 
the data analysis and prediction tools are contributing significantly to current Space Launch 
Initiative turbine airflow test and blade surface pressure prediction efforts. 

Facility Description 

The test series was conducted in MSFC’s cold air flow TTE. 1 The TTE (Fig. 1) is a 
blowdown facility that operates by expanding high-pressure air (420 psig) from one or two-6000 
cubic feet air tanks to either atmospheric or vacuum conditions. (The recent addition of an ejector 
system to the facility extended the facility exit condition capability to vacuum conditions.) Air 
flows from the storage tanks through a heater section, quiet trim control valve, and a calibrated 
subsonic mass flow venturi Flow then continues through the test model, backpressure valve, and 
exhausts to atmosphere or vacuum. The facility can accommodate axial flow, radial inflow, and 
radial outflow turbines. 
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This equipment can deliver up to 220 psia air for run times from 30 seconds to over one 
hour, depending on inlet pressure and mass flow rate. The heater allows a blowdown-controlled 
temperature between 530° R and 830° R. The facility has manual set point closed-loop control of 
the model inlet total pressure, inlet total temperature, shaft rotational speed, and pressure ratio. In 
addition to these control parameters, the facility can accurately measure mass flow rate, torque, and 
horsepower. The associated low-speed data acquisition system is capable of measuring 512 
pressures, 120 temperatures, and several model health-monitoring variables. The high-speed data 
acquisition system will be discussed in a subsequent section. 

Model Description 

The model tested, named the HPFTP Turbine Test Article (TTA), was a full-scale model of 
the Rocketdyne HPFTP turbine with rough rotor blades in the baseline configuration (Fig. 2). The 
model had been tested previously in this configuration, and the configuration along with the 
performance test results are documented in reference 2. Unlike the previous testing, the model 
bullnose was modified for this test to accommodate a slip ring unit for the unsteady pressure 
measurements. 3 The effect of this modification on turbine performance will be addressed when the 
test results are presented. 

Instrumentation 

The TTA contained steady performance instrumentation as well as miniature pressure 
transducers on the blades for fluctuating pressure measurements. The steady state 
instrumentation served two purposes: performance evaluation and model health monitoring. This 
instrumentation included pressures (total and static), temperatures, flow angles, shaft speed 
pickups, and accelerometers. An overview of the steady state model instrumentation is given in 
Table 1. Reference 4 contains more details on the model instrumentation. 
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Instrumented rings defined the turbine inlet and exit planes. Each of the rings 
accommodated a total of eight rakes and two probes with radial actuators. Each rake contained four 
probes (total pressure or Icmperature) at the turbine inlet and five probes (total pressure or 
temperature) at the turbine exit. These rakes could be manually adjusted for yaw angle. The probes 
used with the radial actuatoi s were three-hole cobra probes that were calibrated to obtain yaw angle, 
total pressure, static pressute, and total temperature. These cobra probes were used in the "auto- 
nulling" mode. The circumferential traverse actuators were not used on the inlet and exit rings. 
Both rings were "locked" in ihe 0° position. 

The model included numerous static pressure measurements along the turbine inner 
diameter (ID) and outer diameter (OD) flowpath walls. Static pressure instrumentation was also 
present on the stator vanes along the suction and pressure surfaces at 10%, 50%, and 90% span. 

The 1 st stage turbine blades were instrumented with a total of 24 semiconductor type 
miniature fluctuating pressure transducers manufactured by Kulite. 3 The installed frequency 
response of the transducers was 100 kHz and provided ample bandwidth for the experiment since 
the 1 st stationary vane passage frequency was approximately 4800 Hz. Fluctuating pressure 
transducer footprints were approximately 0.5 mm X 0.5 mm with the sensing diaphragm flush 
with the blade surface (Fig. 3). The sensors were distributed over seven turbine blades at various 
span and chord locations. Table 2 describes the sensor locations in detail. 

Sensor wiring was touted down each blade across the disk to a wire carrier where a pin 
connection was made to a slip ring. 3 A 100 channel high-speed slip ring unit (SRU) routed the 
sensor output to the downsiieam signal conditioning and data acquisition equipment. The model 
was modified to accommodate the SRU and a shaft encoder. Since the model bullnose was 
modified to use the slip ring unit, the total pressure and total temperature rakes on the turbine inlet 


6 



rotating ring were also modified from the previous performance testing. A spacer was used with 
these rakes to lift one of the sensors out of the flow; therefore, each inlet rake only had four probes 
radially in the flow in place of the five used on the previous test. 

High Frequency Data Acquisition 

Banks of remotely controlled Pacific Instrument amplifiers were used to amplify the 
millivolt scale outputs of the on-blade pressure signatures. The unsteady data acquisition system 
provided real-time display and acquisition of all on-blade pressure channels as well as the shaft 
position encoder channels (both once-per-revolution and 500 pulse-per-revolution signatures). 
The 32-channel Computer \ided Dynamic Data Monitoring and Analysis System (CADDMAS) 
was developed as a cooperative research effort with Vanderbilt University, Arnold Engineering 
and Development Center, and Marshall Space Flight Center. 5 6 The CADDMAS is a parallel 
processor based on digital signal processors, analog/digital front-end processors, and standard 
personal computers. Using a parallel processing approach, the system achieves supercomputer 
performance in an interactive, high data-bandwidth environment. 

The real-time capability of the CADDMAS to both display and acquire all of the on-blade 
measurements proved invaluable throughout the test series. Instantaneous waveform display 
aided pre-test sensor calibiations and identified errors in signal conditioning setup prior to test 
runs. With the system, turbine fluctuating pressures were sampled at 85 kHz over acquisition 
sessions lasting on the order of minutes. 

Transducer sensitivities of the installed pressure sensors were provided with delivery of 
the instrumented disk. To ensure transducer calibration accuracy, complete end-to-end (i.e. from 
sensor diaphragm through digital representation) calibrations were performed at least twice each 
test day. Care had to be taken to guarantee that the instrumented blades had reached thermal 
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stability during static step pressurizations before calibration voltages were obtained. A slight 
shift in the bias sensitivities in several of the channels over the duration of the testing was noted. 
Dunn and Haldeman noted similar behavior in their characterization of the SSME fuel turbine. 7 
The researchers attributed most of their long-term sensitivity drift to loss of protective RTV 
coating on the sensor diaphragms. The MSFC turbine model sensors utilized a similar RTV 
protection layer. Future blade surface pressure mapping efforts will utilize a novel temperature 
compensation technique to reduce surface mount pressure sensor error. 8 An order of magnitude 
reduction in temperature induced error is expected through use of the “sense-resistor based” 
compensation technique. 

Test Conditions 

Testing was conducted at a total of 17 set points. First, the turbine’s aerodynamic design 
point from the previous performance test was repeated. 2 This set point is referred to as the “old 
design point” (ODP) or baseline. It corresponds to the SSME 104% rated power level (RPL) 
based on the Rocketdyne engine power balance model of the late 1980’s. Second, the turbine 
conditions were set to match the operating condition in reference 7. This set point is referred to 
as the “Calspan set point” (CSP). Two sets of off-design conditions were then run. First, the 
turbine was set to 65, 70, 80, 90, 100, 104, and 109% RPL based on the updated 1997 version of 
the Rocketdyne engine power balance model. Then, the turbine was run over a range of 
conditions to change the incidence angle on the 1 st stage rotor blades. Incidence angle set points 
included 0, ±5, ±15, and i 25 degrees. Finally, a high turbine pressure ratio set point was run. 
The set point parameters for the test were the turbine inlet total pressure, inlet total temperature, 
speed, and pressure ratio. Table 3 provides a summary of test set point parameters of interest. 
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Test Results and Discussion 


Steady State Performance Results 

Table 4 gives a summary of the steady test results for the old design point (ODP) and the 
Calspan set point (CSP). Data from three test runs done at different times during testing were 
combined to obtain the final numbers for each set point for the current test. Each test run 
consisted of 10 frames of steady data. A sufficient time interval was allowed between each 
frame of data so that the frames could be considered independent results. The inlet and exit total 
pressures and temperatures given in the table are averages of the rake measurements. The turbine 
inlet and exit static pressures are averages of the inner and outer diameter wall static pressure 
measurements at these planes. 

The turbine’s velocity ratio, flow parameter, and efficiency are given in Table 4 to define 
the overall performance. To achieve the proper units, the velocity ratio equation was 
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The turbine flow parameter was defined as follows: 


FP = 
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Note that this is an “engineering” definition for flow parameter that is not truly nondimensional. 

The thermodynamic method of determining turbine efficiency was used. 9 10 The 
temperature drop across the turbine was measured to determine the actual enthalpy change. For 
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this “cold” turbine testing where the temperature was relatively low, an ideal gas was assumed 
and y and C p were considered constant. 
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Turbine steady-state performance across both design and off-design conditions is shown 
in Fig. 4. This figure shows the thermodynamic turbine efficiency versus velocity ratio. 

A detailed posttest uncertainty analysis was completed for the current test using the 
methodology in reference 11. Estimates for the systematic and random components of the 
uncertainties of the measured variables are given in Table 5. The random component estimates 
were obtained directly from multiple test results based on a large sample assumption, and they 
include set point repeatability. 

Data at the old design point test condition was compared to data obtained from previous 
testing to evaluate the effect of the unsteady instrumentation on the turbine’s performance Static 
pressure drops through the turbine as well as overall performance parameters were studied. 

Figure 5 gives the static pressure drops through the turbine for both the current test and 
the previous test. 2 The y-a\is in Fig. 5 is the average of the inner and outer wall static pressure 
measurements at each plane normalized by the static pressure at the turbine inlet. The \-axis 
represents each axial station: 1 is the turbine inlet, 2 is the l sl stator inlet, 3 is the 1 st stator exit, 4 
is the 2 nd stator inlet, 5 is the 2 nd stator exit, 6 is the turbine exit, 7 is the EGV inlet, and 8 is the 
EGV exit. Figure 5 shows that there was no measurable difference in the static pressure drop 
through the turbine between the current and previous tests. 
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To compare overall performance, several parameters were studied. These parameters 
included facility set points as well as calculated performance parameters. The facility set points 
were Poi, T 0 i, N, and Pr; therefore, velocity ratio was a set point. The calculated performance 
parameters that were used to compare the two data sets included the flow parameter and 
thermodynamic efficiency. The data from both the current test and the previous test is given in 
Table 4. Estimates for the random component of the uncertainty of the calculated values are 
given in Table 6 for the current test. The random component estimates were obtained directly 
from multiple test results based on a large sample assumption, and they include set point 
repeatability. Values for the previous test were similar. 

Table 4 shows that t he pressure ratios and the velocity ratios were the same for both tests. 
The flow parameter increased by 2.4% for the current test. The subsonic venturi used to measure 
the mass flow rate in the TTE was recalibrated between tests. The application of the new 
calibration data caused a shift in the measured mass flow rate and the corresponding flow 
parameter. 

Thermodynamic turbine efficiency was compared between the two tests. The calculated 
thermodynamic efficiency was 85.7% for the previous test and 85.3% for the current test. This 
gives Ar(=0.4%. The methodology in reference 12 was used to evaluate the difference in 
efficiency obtained from the two tests. The temperature measurements for both tests were made 
with the same thermocouples, wiring, hook-up, and data acquisition system. Similarly, the 
pressure measurements were made with the same instruments, hook-up, and data acquisition 
system. No factors were identified to change the systematic uncertainties of the measurements 
used to calculate the thermodynamic efficiency between the two tests. The random component of 
the uncertainty in the thermodynamic efficiency, including test-to-test variations, was 0.0014 or 
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0.14%. Therefore, the uncertainty in the difference in efficiency between the two tests was 
calculated to be 0.004 or 0.4%. 12 Since Ar| ± U Ar 1 includes 0, the difference is insignificant. 12 

In summary, the dau comparison showed no significant changes in turbine performance 
between the two tests. This leads to the conclusion that the installation of the unsteady 
instrumentation and the modification to the model bullnose did not have a measurable effect on 
the turbine’s performance. 

Unsteady Performance Results 

Several of the turbine blade surface mounted pressures failed prior to the start of testing. 
Of the original 24 transducers, only 18 were operational for the first instrumented rotor test with 
three more lost during testing. In interpreting the high frequency turbine blade surface pressure 
data, both standard power spectral density (PSD) and synchronous time averaging (STA) 
processing were employed.' ’ 141516 The high-speed data sample rate across all testing was set at 
85334 Hz which was sufficient to capture the sixth harmonic of the 1 st stationary vane passing 
component. The PSD block size for the data archival was set to 8192 yielding a resolution 
bandwidth of 10.4 Hz, or approximately 0.23% of the I st vane passing rate at baseline set point 
conditions. A slight periodic rotor speed variation (measured in revolutions per minute, RPM) 
was noted in testing, but did not affect the blade surface pressure characterization effort. 
Figure 6 displays the periodic variation in frequency of the 1 st vane passing rate as well as its first 
harmonic for the baseline set point at 6982 RPM. The variation was approximately 18-RPM 
peak-to-peak at 0.19 Hz. STA processing provided pressure waveforms free of random and 
nonsynchronous disturbances. The STA process utilizes a trigger-pulse (supplied by the shaft 
key phaser) to synchronize the averaging process with shaft rotation. The resulting waveform 
over one full revolution contains information that is purely synchronous related. The signal 
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analysis technique is equivalent to a comb filter with center frequencies at the synchronous rate 
and subsequent harmonics Both the maximum peak-to-peak variation and mean of these 
pressure waveforms were used as primary blade surface pressure assessment attributes. 

All displayed unsteady results have been normalized by model inlet total pressure. 
Negative chord values denote pressure surface placement while positive denote suction surface. 
In terms of percentages, the leading edge of the turbine blade corresponds to a 0% chord, while 
100% chord denotes trailing edge. Averaging frame length for all unsteady parameters was 
25 seconds, which guaranteed at least 980 shaft revolutions per estimate at the lowest test set 
point speed of 2344 RPM. 

A sliding-average process was applied to STA waveforms to eliminate what appeared to 
be slip-ring related noise. The first equivalent low-pass filter zero-crossing for the averaging was 
approximately 195N. With a 1 st stage vane rotor count of 41, the pass band contained the first 
three harmonics of 1 st stationary vane passing. 

Figure 7 shows the normalized average surface pressure distribution along the chord at 
blade mid-span for baseline conditions. The variation across test repeats reflects sensor bias 
drifting which could not be removed from the acquired data. Implementation of the calibration 
method of Ainsworth 8 is expected to significantly improve the thermally driven uncertainty in 
future turbine blade pressure mapping efforts. Figure 8 shows the normalized maximum 
peak-peak variation (or unsteady envelope) of surface pressures also at mid-span across chord for 
the baseline condition. The unsteady envelope was very repeatable across the baseline runs as 
shown across pressure locations. 

Figure 9 displays the chordwise unsteady pressure distribution for 50% span across 1997 
SSME Power Balance Rated Power Levels (RPLs); i.e., from 65% to 109% RPL. Interestingly, 
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the fluctuating pressure variation across power level is only slightly greater than the test-to-test 
variation recorded at the baseline set point (Figure 8). This suggests minimal high-speed fuel 
turbine blade load variation across the SSME flight envelope. Figure 10 shows the 50% span 
normalized unsteady pressure distribution across incidence variation from -25 degrees to 
+25 degrees. Peak leading edge fluctuating pressures for the entire experimental series occurred 
at the +25 degree incidence condition. 

Detailed time-resolved pressure waveforms will be shown and discussed in the 
comparison of experimental data to CFD prediction. 

Unsteady CFD 

Computational Procedure 

The computational tool used to simulate the flow in the turbine is known as CORSAIR. 
The governing equations solved in CORSAIR are the time dependent, three-dimensional 
Reynolds-averaged Navier-Stokes equations. The algorithm used in the computational procedure 
consists of a time-marching, implicit, finite-difference scheme. The procedure is third-order 
spatially accurate and second-order temporally accurate. The inviscid fluxes are discretized 
according to the upwind scheme developed by Roe. 17 The viscous fluxes are calculated using 
standard central differences. An approximate-factorization technique was used to compute the 
time rate changes in the conserved variables. The turbulent viscosity is calculated using a highly 
modified two-layer Baldwni-Lomax algebraic turbulence model. The numerical analysis uses 
Message Passing Interface (MPI) and OpenMP to reduce the computation time for large-scale 
three-dimensional simulations. 

The Navier-Stokes analysis uses O- and H-type zonal grids to discretize the flow field and 
facilitate relative motion of the rotating components. The O-grids are body-fitted to the surfaces 
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of the airfoils and generated using an elliptic equation solution procedure. They are used to 
properly resolve the viscous flow in the blade passages and to easily apply the algebraic 
turbulence model. The algebraically generated H-grids are used to discretize the remainder of the 
flow field. Figure 1 1 shows the computational grid at mid-span for the SSME geometry. 

Boundary Conditions 

For subsonic inlet I low, the total pressure, total temperature, and the circumferential and 
radial flow angles are specilied as a function of the radius. The upstream-propagating Riemann 
invariant is extrapolated from the interior of the computational domain. 

For subsonic outflow the circumferential and radial flow angles, total pressure, and the 
total temperature are extrapolated from the interior of the computational domain. The ratio of the 
exit static pressure to the inlet total pressure is specified at mid-span of the computational exit, 
and the pressure at all other radial locations at the exit is obtained by integrating the equation tor 
radial equilibrium. Periodicity is enforced along the outer boundaries of the H-grids in the 
circumferential direction. 

No-slip boundary conditions are enforced along the solid surfaces. It was assumed that 
the normal derivative of the pressure was zero at solid wall surfaces. All solid surfaces were also 
assumed to be adiabatic. The flow variables at zonal boundaries were explicitly updated each 
time step by interpolating values from the adjacent grids. 

Geometry. Grids, and Simulation Characteristics 

The two-stage turbine geometry contained 41 first-stage vanes, 63 first-stage rotors, 39 
second-stage vanes, and 59 second-stage rotors. A 40-60-40-60 blade count approximation was 
made allowing 2-3-2-3 simulations. To keep the pitch-to-chord ratio constant, the first-stage 
vane airfoils were scaled by factor of 41/40, the first-stage rotor blades were scaled by a factor of 
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63/60, the second-stage vane airfoils were scaled by a factor or 39/40, and the second-stage rotor 
blades were scaled by a factor of 59/60. The rotor tip clearances were set to 1.4% of the rotor 
height. A total of 3,556,785 grid points were used. The average value of y + , the non- 
dimensional distance of the first grid point above the wall, was approximately 1.5 for the airfoil 
surfaces and 3.5 for the endwalls. 

The simulations were performed on 20 400-MHz processors of an 0rigin2000 located at 
NASA Ames Research Center. The 20 processors consisted of 10 MPI processes and two 
OpenMP threads per MPI process. The simulations were run for 10 characteristics cycles, where 
a characteristic cycle consists of three rotor blades passing by two vane airfoils. The simulations 
were performed at 12,000 iterations per cycle and required approximately 3x10° sec/grid 
point/iteration computation time. 

Results 

The predicted time- averaged flow values and performance parameters for the “baseline 
set point” equivalent turbine inlet conditions are presented in Table 7. The flow values were 
normally recorded midwav between blade rows, except for the first row where the values were 
recorded at the computational inlet (approximately one axial chord upstream of the first-stage 
vane leading edge) and the- last row where the values were recorded at the computational exit 
(approximately one axial chord downstream of the second-stage rotor trailing edge). The overall 
total -to-total efficiency of die turbine was 0.888 based on the change in the total temperature and 
0.870 based on the change in circumferential velocity. The overall total-to-static efficiency of 
the turbine was 0.843. 

Instantaneous entropy contours at the mid-span section of the turbine for the design 
operating conditions are shown in Fig. 12. Entropy contours are useful in highlighting the 
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convection of airfoil wakes- and regions of high loss. Figure 12 clearly shows the stretching of 
the vane wakes as they pass through the downstream rotors. The stretching of the wakes is 
caused by the differences m the core velocities on the pressure and suction sides of the rotor 
passages. The airfoil wakes can persist through several blade rows, leading to an accumulation 
of wakes in the downstream rows. In a subsonic turbine, such as the one under consideration in 
this study, the airfoil wakes tend to drive the dominant unsteady frequencies in adjacent blade 
rows. The accumulation of wakes in the downstream rows can also lead to the presence of 
unsteadiness at harmonics of the fundamental frequencies. 

Comparisons with Experiment 

For steady state turbine performance, CFD predicted an overall total-to-total efficiency of 
0.888 for the baseline set point. The corresponding experimentally derived efficiency for the 
current test (Table 4) was 0.853. It is not unusual for CFD simulations to yield a higher overall 
efficiency; the CFD model does not contain all the loss mechanisms present in the experimental 
apparatus. The CFD mode l includes only the viscous losses due to the airfoil boundary layers, 
the endwall boundary layers (and associated secondary flows), and the tip clearance flows. 

Figure 13 contains predicted versus experimental normalized time-averaged surface 
pressure distributions (baseline set point) at 50% span of the first-stage rotor while Figure 14 
contains the unsteady distributions. All values have been normalized by inlet total pressure. The 
CFD results exhibit good agreement with both the experimental time-averaged values and the 
amplitude of the unsteadiness. 

Figures 15 and 16 compare the instantaneous pressure waveforms for arbitrary pressure 
and suction surface locations on the first-stage rotor, respectively. Instantaneous pressure is 
plotted versus rotation index where an index of 0.1 represents 10% of a shaft revolution. In 
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developing the overlays, the experimental and computed waveforms were arbitrarily aligned to 
match phase. As can be seen in the comparisons, the CFD prediction not only matches 
fluctuating pressure amplitude, but it also predicts the shape of the waveform reasonably well. 
Both the experimental and predicted results indicate that the predominant unsteadiness is 
associated with the periodic passing (as seen from the rotor) of the upstream vane wakes. The 
amplitude of the unsteadiness is greater on the suction surface of the rotor due to its increased 
exposure to the vane wakes and potential field. 

In addition to the baseline set point, a CFD simulation was also performed for an 
off-design turbine condition corresponding to +15 degrees incidence on the first-stage rotor. 
Figure 17 contains the predicted and experimental normalized time-averaged surface pressure 
distributions at 50% span, and Figure 18 contains the corresponding unsteady pressure envelopes. 
While the predicted time averaged surface pressure distribution compares well with the 
experimental data, the suction-side leading edge unsteady distributions exhibit some 
discrepancies. In particular, the predicted suction side distribution is shifted more towards the 
leading edge than is indicated by the experimental data. The main sources of these discrepancies 
are (1) The experiments were performed with rough blade surfaces while the CFD simulations 
assumed smooth blades. This difference has no effect on the set point results but does affect the 
incidence angle at which the suction surface boundary layer will begin to separate. (2) The 
turbulence model in the CFD code, although well anchored for turbomachinery simulations, will 
begin to break down at far off-design flow conditions. 

Figures 19 and 20 compare the instantaneous pressure waveforms for the same arbitrary 
pressure and suction surface locations shown for the baseline configuration. The predicted 
pressure surface trace displays excellent agreement with the experimental data. Both the 
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predicted results and experimental data indicate that the predominant unsteadiness is associated 
the vane-passing frequency. The experimental trace on the suction surface shows the 
predominant unsteadiness occurs at the vane-passing frequency, while the predicted results also 
contain higher-frequency content. This additional frequency content is consistent with 
intermittent breakdown of the suction-surface boundary layer in the simulations due to the 
assumption of smooth blades and limitations of the turbulence model. 

Summary and Conclusions 

A technology program to increase the understanding of unsteady environments for rocket 
engine turbines was successfully completed. Steady measurements used to evaluate turbine 
performance as well as unsteady pressure measurements on the turbine rotor blades were 
obtained. 

The steady measurements were used to define the turbine’s performance at each set point. 
Additionally, the existence of a comprehensive steady data set on the turbine tested allowed a 
unique opportunity to evaluate the effect of the installation of the unsteady instrumentation on 
the turbine’s performance. Data comparisons between the current test and a previous test of the 
same configuration with no unsteady instrumentation showed that the unsteady instrumentation 
and resulting hardware mollifications did not affect the performance of the turbine. This data 
also demonstrated excel leni facility repeatability, both run-to-run and test-to-test. 

With the unsteady data, success in the acquisition and evaluation of complex high-speed 
turbine on-blade surface pressure environments was shown. The benefit of advanced high-speed 
digital signal acquisition arid processing was demonstrated in the accurate mapping of the blade 
unsteady pressures. The time-averaged and unsteady CFD results showed very good agreement 
with the experimental data at the design operating conditions and fair-to-good agreement at an 
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operating condition corresponding to +15 degrees of incidence on the first-stage rotor. The 
high-fidelity flow features characterized in the successful effort will undoubtedly aid the 
advancement of CFD prediction tools. Lessons learned will be invaluable in the planning and 
implementation of similar blade mapping efforts for advanced turbines for the next generation 
and future generation reusable launch systems. 
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Fig. 2. Turbine Test Article 
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Fig. 3. Instrumented Blade (50% Span Suction Side) 
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Fig. 4. Thermodynamic Efficiency vs. Velocity Ratio 
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Fig. 5. Normalized Static Pressure Drop Through Turbine 
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Fig. 6. Shaft Speed Variation 
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Fig. 8. Unsteady Surface Pressure Envelope at 50% Span for Baseline Conditions 
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Fig. 11. Mid-span Computational Grid 
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Fig. 12. Mid-span Instantaneous Entropy Contours at the Design Flow Conditions 
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Fig. 13. Time-averaged Surface Pressure at 50% Span versus CFD Prediction — Design Flow 
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Fig. 14. Unsteady Envelope Pressure at 50% Span versus CFD Prediction — Design Flow 
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Fig. 17. Average Surface Pressure at 50% Span versus CFD Prediction at Off-Design 
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Fig. 18. Unsteady Envelope Pressure at 50% Span versus CFD Prediction at Off-Design 
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Table 1. Steady State Instrumentation Overview 


Turbine Inlet and Exit: 

4 total pressure rakes (4 probes each at inlet and 5 probes each at exit). 

4 total temperature rakes (4 probes each at inlet and 5 probes each at exit). 

2 auto-nulling cobra probes with radial actuators. 

Automatic circumferential traverse. 

Turbine: 

8 inner and 8 outer wall static pressures at 6 axial planes. 

Stator surface static pressures — 6 on pressure side and 8 on suction side at 10% span, 50% span, 
and 90% span on both stages. 

14 stator outer shroud and 14 stator inner shroud static pressures on both stages. 

Disk cavity static pressures. 

Disk cavity total temperatures. 

Exit Guide Vanes: 

12 inner and 12 outer wall static pressures at 2 axial planes. 

4 total pressure measurements on 6 vanes. 

4 total temperature measurements on 6 vanes. 

Miscellaneous: 

2 speed pick-ups. 

Accelerometers — 2 horizontal, 2 vertical. 

Contoured blank plugs for all bosses. 

Health monitoring instrumentation. 
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Table 2. On-Blade Instrumentation Locations 


Blade 

# 

% Span 

% Axial Chord 

% Wetted Length 

Suction/ 

Pressure 

21 


14 

12.5 

P 

21 

90 

58 

47.6 

P 

21 

90 

75 

65.5 

P 

22 

90 

7 

8.5 

S 

22 

90 

22 

19.1 

S 

22 

90 

88 

81.3 

s 

42 

10 

13 


p 

42 

10 

59 

49.6 

p 

42 

10 

77 


p 

43 

10 

5 

8.3 

s 

43 

10 

18 


s 

43 

10 


51.0 

s 

43 

10 

83 

74.7 

s 

63 

50 


9.4 

p 

63 

50 

23 

20.3 

p 

63 

50 

52 

44.1 

p 

63 

50 

82 

75.0 

p 

1 

50 

5 

7.3 

s 

2 

50 

9 

11.5 

s 

1 

50 

19 

19.0 

s 

2 

50 

25 

24.0 

s 

1 

50 

61 

50.1 

s 

2 

50 

72 

61.9 

s 

1 

50 

89 

83.0 

s 
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Table 3. SSME HPFT Test Matrix 


Run/Rerun 

Poi (psia) 

Toi (°R) 

N (RPM) 

Model Pr t .t 

Description 

8/0 

100 

550 

6982 

1.47 

Old Design Point 

8/1 

100 

HBLkiiW 

■jgEEM 

1.47 

Old Design Point 

8/2 

Fob 

550 

| 

1.47 

Old Design Point 


50 

550 

6747 

1.61 

Calspan Set Point 

9/1 

■SB 

550 

6747 

1.61 

Calspan Set Point 

8/3 

— 

550 

6982 

1.47 

Old Design Point 

11/0 

100 

550 

6791 

1.51 

97 PB_104%RPL 

12/0 

100 

550 

6979 

1.52 

97 PB_109%RPL 

13/0 

— 

550 


1.38 

1 ■ 1 1 II 111 

14/0 


550 



97 PB_70%RPL 

15/0 

100 

550 

5970 

1.43 


16/0 

100 

550 

6306 

1.46 

97 PB_90%RPL 

17/0 

100 

550 

6649 

1.49 

97 PB_100%RPL 

18/0 

100 

550 

2344 

1.51 


19/0 

100 

550 


1.51 



100 

550 

4895 

1.51 

+15 deg incidence 


100 

550 

8032 

1.51 

-15 deg incidence 


100 

550 

6142 

1.51 

+5 deg incidence 


100 

550 

7182 

1.51 

-5 deg incidence 


100 

550 

6728 

1.51 

0 deg incidence 

— 

100 

580 

7951 

1.84 

High Pr Point 


M 

580 

7951 

1.84 


9/2 

SB 

550 


1.61 


8/4 

100 

550 

6982 

1.47 

Old Design Point 

9/3 

50 

550 

6747 

1.61 

Calspan Set Point 
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Table 4. Steady State Performance Results 



Previous Test 

Current Test 


ODP 

ODP 

CSP 

Facility Measurements 



99.90 

50.26 

To Fac 


548.00 

548.16 

N 

7004.45 

7005.98 

6763.11 

W 

13.98 

14.37 

7.73 

Model Conditions 

Poi 

99.47 

99.89 

50.24 

T 01 

547.64 

545.71 

545.90 

P. 

99.30 

99.65 


Po2 

67.67 

67.94 

31.29 

To2 

498.65 

497.17 

488.53 

P 2 

66.34 

66.66 

30.43 

Pressure Ratios 

t-t 

1.47 

1.47 

1.61 

t-s 

1.50 

1.50 

1.65 

Overall Performance 

U/C 0 

0.37 

0.37 

0.33 

FP 

3.29 

3.36 

3.59 

Tlth 

0.857 

0.853 

0.831 
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Table 5. Uncertainty Estimates of Measured Variables 



ODP Random 

CSP Random 

Systematic 

Po Fac 


fflSSElBMKKM 



0.10° R 

0.27° R 

1.0 °R 

N 

1 .3 RPM 

1.3 RPM 

1.0 RPM 

W 

0.0037 lbjs 

0.010 lbjs 

1% reading 

Po. 


Q3SBIHI 


To, 

0.11° R 

0.26° R 

o 

0 

JO 

P. 

0.034 psi 

EEnMi 


P 02 

0.057 psi 



To2 

0.079° R 

0.16° R 

Oi 

o 

o 

P 2 
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Table 6. Random Component Uncertainty Estimates of Calculated Values 



ODP Random 

CSP Random 

Pr t -t 

0.0011 

0.0023 

Pr ,-s 

0.0012 

0.0024 

U/C 0 

0.0004 

0.0005 

FP 

0.0011 

0.0078 

Tlth 

0.0009 

0.0008 
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Table 7. Predicted Time-Averaged Flow Values and Performance Parameters at Design 

Conditions 



Vane-1 

Rotor- 1 

Vane-2 

Rotor-2 

M, (abs) 

0.138 

0.432 

0.178 

0.430 

M, (rel) 

0.138 

0.203 

0.178 

0.205 

M2 (abs) 

0.432 

0.178 

0.430 

0.179 

M2 (rel) 

0.432 


0.430 

0.355 

Pi (abs) psia 

98.68 

86.87 

80.44 

81.29 

P2(abs) psia 


80.44 

81.29 

66.67 

P01 (abs) psia 


98.88 

82.28 

81.29 

P01 (rel) psia 

100.0 

89.48 

82.28 

73.69 

P02 (abs) psia 

98.88 

82.28 

81.29 

68.19 

P 02 (rel) psia 

98.88 

88.40 

81.29 

72.79 

T 0 i (abs) R 

546.2 

546.2 

520.0 


T 0I (rel) R 

546.2 

530.8 

520.0 


T02 (abs) R 

546.2 

520.0 

520.0 

495.9 

T 02 (rel) R 

546.2 

530.8 

520.0 

505.3 

a, (deg) 

0.000 

-68.72 

14.56 

-67.13 

Pi (deg) 

0.000 

- 37.02 

14.56 

-32.00 

a 2 (deg) 


14.56 

-67.13 

6.70 

P2 (deg) 

wmmxnmm 

62.86 

-67.13 

60.16 

W (lbm/sec) 

15.2 

15.2 

15.1 

15.1 

Power (hp) 

-- 

136 

— 

123 

fll-t (based on AV) 

— 

0.874 


0.867 

(based on ATt) 

— 

0.884 

— 

0.887 

Tlt-s 

— 


— 

0.795 

Zweiffel 

0.508 

1.145 

0.625 

1.188 
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